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ABSTRACT 

Positions of High-Mass X-ray Binaries are often known precisely enough to unam- 
biguously identify the optical component, and a number of those stars are monitored 
by the OGLE and MACHO collaborations. The light curves of two such candidates 
are examined for evidence of Be star behaviour and for periodicity. One of the stars 
exhibits two periods of 6.833 and 15.586 hours, much shorter and more stable than 
periods of Be/X-Ray Binaries that are attributed to the Be star's disk, but consis- 
tent with short-term Be variability attributed to pulsations. The multiperiodicity is 
quantified with Fourier techniques and examined for phase stability; a combination of 
radial and non-radial pulsations is discussed. 

Key words: stars: oscillations - stars: emission line. Be - X-rays: binaries - methods: 
data analysis 



1 INTRODUCTION 

Due to heavy coverage by recent X-Ray space missions, the 
number of known High-Mass X-ray binaries (HMXB) in the 
Small Magellanic Cloud (SMC) is rapid ly expanding, and 
is now in the sixties (Ha berl et al.ll2nn4al ). The more recent 
space missions, Chandra and XMM-Newton, are able to de- 
liver positions accurate to a few arcseconds, opening the 
door to unambiguous optical identifications. 

Meanwhile, in the optical community, dense starfields 
are being monitored every few nights for variability with the 
goal of detecting massive compact objects via gravitational 
microlensing. Many other areas of astronomy likewise bene- 
fit from well-sampled light curves in the optical and near-IR. 
The MACHO project ran for seven years, until 1999, and has 
recently put its data online as a service to the community. 
The OGLE project is ongoing, is currently in its third stage, 
and is expanding the area covered with each successive re- 
alization. 

Careful analysis of a number of these optical counter- 
parts to X-ray sources, which are frequently Be stars, is 
yielding: interest i ng mu ltiperiodic behaviour. For instance, 
ISchmidtke~ give four examples of stars which 

have an orbital period, judged by outbursts, of order one 
hundred days as well as quasi-periodic variations of order 
10 days. A region on a disk excited by the passing compact 
obj ect may cause t his latt er variation as it orbits the star. 

iHaberl et al.1 (|2nn4bl) identified two new X-ray pul- 
sars with bright st ars in the OGLE catalog of the SMC 
(jUdalski et al.lll998i) . I examine these stars for periodic vari- 
ations, finding no result for the star associated with the 499.2 
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s pulsar XMMU J005455.4-724512 ("star A") but finding 
two periods that are both less than one day for the star as- 
sociated with the 701.6 s pulsar XMMU J005517.9-723853 
("star B"). In section 2 I explain how the light curves were 
constructed. A search for significant periods is covered by 
section 3. Section 4 deals with the evolution of those peri- 
ods over the time-scale of the observations. I close with a 
discussion of the possible physical mechanisms causing this 
multiperiodicity in section 5. 



2 LIGHT CURVES 

I obtained OGLE data for both stars and for nearby cali- 
bration stars from Andrzej Udalski. This paper is the first 
time the OGLE light curves for these two objects have 
been published, since their amplitudes of variability are too 
small for merabersh ip in the catalog of SMC variable stars 
feeb run et ani2nnih . They include both OGLE H and the 
ongoing OGLE HI data, with the break at MJD 52000. At 
the present time OGLE HI is not carefully calibrated, so 
I selected stars within 10 arcseconds of the stars of inter- 
est that were brig:hter than magnitude 20 in the OGLE II 
catalog (' Udalski et al.l ll998^ to make a crude calibration. I 
compared the difference in average magnitude of each star 
across the OGLE II/OGLE HI break, then determined an 
offset to apply to the break by inverse variance weighting 
those individual offsets. These total offsets have errors of 15 
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Figure 1. These are the V and E MACHO and / OGLE hght 
curves for star A, which is associated with XMMU J005455.4- 
724512. Errors are omitted for clarity. 
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Figure 2. Same as above, but for star B, which is associated with 
XMMU J005517.9-723853. 



and 32 mmag for star A and B respectively, which is about 
the same as both the day-to-day and year-to-year variabil- 
ity, so more careful calibration will not significantly affect 
these results. 

The MACHO light curves are publicly available via an 
interactive browser^. I derived standard Fand R magnitudes 
from the instr umental m agnitudes of MACHO 's blue and red 
bands as Alco ck et al.l (p.999) instructs. 

The resultant light curves are plotted in Figs. 1 and 
2. For the MACHO data, some anomalous points are not 
within the plot range. These figures show that the MACHO 
and OGLE times of observation overlap between MJD 50650 
and MJD 51546.5, almost three seasons of data. The rela- 
tive strength of the periodicities during this time will be 
compared and discussed in terms of a physical interpreta- 
tion. 

For the following analysis I used the MACHO blue and 
red bands, not the standard bands. The transformation re- 
quired good data points for both bands simultaneously; by 
using the instrumental values I was able to retain the obser- 
vations for which only one band had good data. 



3 PERIOD ANALYSIS 

As a first attempt to find period s, I made a L omb peri- 
odogram using the modification of lScargl3 (|l982l) . 

Fig. 1 shows changes in the brightness of star A with 
a time-scale of the length of the data, but whether it is 
periodic obviously ca nnot be de termined. Using the false- 
alarm probability from lHorne B aliunas (1986) to measure 
significance, none of the datasets yielded significant periods 
besides the power at a few hundr e d days, which can be seen 
in the light curve. H aberl et al.l (|2004bf) proposed it as a 
Be star because it is within the error circle of an H-alpha 
emitting object (|Meyssonnier (fe Azzopardjll^Sj) • The light 
curve wanders on many time-scales, which is characteristic 
of Be variability, but it has a smaller amplitude than most Be 
stars. It is possible that the H-alpha emission stopp ed after 
the observations of Mev ssonnier Azzopardil fl993), which 
were done in August of 1982, but before MACHO started 
taking data, so the photometric variations may represent the 
star and disk in a quiescent state. 

For star B, the OGLE data showed a dominant period 
of 0.6494 it 0.0001 days (1.5398 ± 0.0003 cycles/day, here- 
after /i). Next I analyzed the MACHO data, which has 
peaks clustered around 0.2847 days (3.5119 ± 0.0010 c/d), 
hereafter /2. On closer inspection it became clear that both 
data sets had the other periodicity too, but it was lost in 
the aliases of the stronger period. 

The spectra of evenly sampled data exhibit aliases at 
the significant frequency plus and minus multiples of the 
sampling frequency, and one can only reliably determine pe- 
riods up to the Nyquist frequency, half the sampling fre- 
quency, keeping in mind that aliases from higher frequencies 
can appear. For unevenly sampled data, aliasing does not 
render the periodicity ambiguous. In that case, the true pe- 
riod can usually be gleaned from the peak of the spectrum 
since there is no solid sampling frequency. Of course, for 
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Figure 3. Star B's periodogram for MACHO blue (top) and red 
(middle). The dotted line is 1% false-alarm probability and the 
dashed line is 0.1% false- alarm probability. The spectral window 
(bottom) is very nearly the same for both. 
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Figure 5. Top to bottom: MACHO blue and red, and OGLE / 
band CLEANed spectra. 
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ground-based photometry, observations can only be taken at 
night, and this sampling introduces secondary peaks in the 
window function at multiples of one day. If the telescope is 
in a location from which the star cannot be observed year- 
round, another small year~^ = 0.002738 c/d peak in the 
window function occurs. Finally, the total length of the ob- 
servations sets the width of a spectral feature, since two 
trial periods whose associated frequencies are closer than 
that will not go out of phase during the observation. 

Fig. 3 shows the periodogram and spectral window for 
the MACHO data of star B. Plotted is the modulus squared 
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Figure 6. The peak-to-peak amplitude of the periodic variation 
near 1.54 c/d, plotted by year for all data sets. MACHO modu- 
lation in this region was not statistically significant during 1993 
and 1995, so these years are omitted. Errors were determined by 
finding the variance in the peak-to-peak amplitude of an ensem- 
ble of light curves created by varying each datum by a Gaussian 
of standard deviation equal to its error. MACHO blue data are 
empty circles, MACHO red data are filled circles, and OGLE data 
are crosses. 
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Figure 7. Same as above, near 3.51 c/d instead. OGLE modu- 
lation in this region was not statistically significant during 1998, 
so this year was omitted. 



of the spectral window; at negative frequencies reside the 
complex conjugate of the positive values. Although bad 
points were removed only from the band in which they were 
bad, both bands share the same time sampling, so the win- 
dow function is nearly identical. As explained, the spectral 
window has its strongest aliases at multiples of 1 c/d, but 
they are weaker than the true value, so naively picking the 
peak as the true value works in this case. Next to the /2— 2 
c/d = l.bl c/d alias lies another significant peak at /i = 1.54 
c/d. Its aliases can also be seen. Fig. 4 is the periodogram 
and spectral window for the OGLE data for star B. Here 
the fi peak is dominant with the /2 peak hiding among its 
aliases. 

The CLEAN algorithm was designed to extract multiple 



Figure 8. /I's frequency by year, judged from the peak of the 
Lomb-Scargle periodogram. MACHO modulation in this region 
was not statistically significant during 1993 and 1995, so these 
years are omitted. MACHO blue data are empty circles, MACHO 
red data are filled circles, and OGLE data are crosses. Typical 
FWHM is 0.004 c/d. 
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Figure 9. Same as above, but for /2. OGLE modulation in this 
region was not statistically significant during 1998, so this year 
was omitted. 



periodicities from uneve nly sampled data de spite its compli- 
cated window function (jRoberts et al.lll987l) . I ran CLEAN 
on all three datasets and discovered both frequencies present 
in all of them (see Fig. 5). For each light curve the periodic- 
ities have different powers; this could be due to differences 
in both bandpass and epoch of observation. In what follows 
I attempt to deconvolve this dual dependence. 

I split the data into year-long chunks and analyzed them 
separately. Due to its position at Las Campanas Observa- 
tory, Chile, the OGLE project is not able to monitor stars 
in the SMC all year round. Most sets of observations began 
around June and went through February; for convenience I 
labeled such seasons with their starting year. For instance, 
2003 in Fig. 6-12 went until February, 2004. I broke the MA- 
CHO data into similar chunks. For each of these chunks, I 
recorded the peak frequency and amplitude in the vicinity 
(within 0.01 c/d) of the two modulations if it was significant 
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Figure 10. OGLE data for star B with power removed between 3.5023 and 3.5223 c/d, then folded at 1.5398 c/d. Each plot is one year 
of data, starting in 1997. The last plot is all seven years of data. Errors are omitted for clarity. 
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at a — 0.01 accord i ng to the false-alarm probability test of 
iHorne (fe BaliimaH U^E^. Fig. 6 and 7 are the two mod- 
ulations' amplitude plotted versus year. For /2, the bluer 
colours show considerably larger modulations, which is true 
for fi to a lesser degree. See the discussion section for a pos- 
sible interpretation. The peak frequencies for each year are 
Fig. 8 and 9. 

It is also instructive to view the folded light curves 
for each year. If one plots the OGLE data folded at /2, 
a faint pattern is discernible by eye. If, however, the data 
is prewhitened with frequencies near /i, the modulation is 
clear. The CLEAN algorithm sequentially identifies and re- 
moves the strongest frequency, so as they were reported, if 
they fell within 0.01 c/d of the dominant /i variation, I re- 
moved the best-fitting sinusoid of that frequency from the 
time series. Fig. 10 and 11 show each year's OGLE data af- 
ter this procedure: they are folded at each frequency after 
removing the other one. 



4 PHASE STABILITY 

These yearly light curves, cleaned of the other period, also 
help probe the period stability of each variation. I found that 
both periods were stable enough so that period stability can 
be described well as a change in phase with respect to a sine 
curve with a fixed frequency. 

Fig. 12 has phase diagrams for each periodicity. They 
were constructed by splitting up the yearly data into 
pieces of about 20 sequential observations. The pieces never 
spanned OGLE's off-season between February and June. 
Next, each piece was folded and a sinusoid was fit to it by 
least squares. Plotted is the phase of that best-fitting sinu- 
soid versus the mean time of observation. A Monte-Carlo 
error estimate for that phase is also plotted: I varied each of 
the data points by a Gaussian of standard deviation equal to 
the quoted error and recalculated the phase, thus building 
up a list of 100 phases. The standard deviation of that list 
is what I adopted as a measure of the error in phase. Large 
error bars mean that offsetting the data had a large effect 
on the best fit, so in these cases it can be understood that 
the sinusoid fit the data only poorly from the start. How- 
ever, if the peak-to-peak amplitude of the fit was less than 
10 mmag, I judged the phase unreliable and did not plot it. 

To give the plot continuity, I constrained the adjacent 
points in time for each band to have phase difference less 
than 0.5 cycles. It should be noted that the phases plot- 
ted may be different by multiples of 1 cycle from reality. 
For instance, in the 1.5398 c/d phase diagram in Fig. 12, 
year 1996 may not be properly matched to year 1995 due to 
sparse data and a quickly changing phase, and in the 3.5119 
c/d diagram, the first data points for each MACHO band 
may be plotted 1.0 cycle too high. 

It is interesting to compare these phase diagrams with 
Figs. 10 and 11. The bottom left box of Fig. 10 corresponds 
to the 2001 data. For 2001 the peak of the periodogram is at 
1.5430 c/d instead of /i = 1.5398c/d (see Fig. 8). This period 
decrease induced an advancing phase, so the folded plot does 
not look as crisp as other years. If we consider the variation 
as simply a sinusoid with a discretely changing period, the 
change of about 0.003 c/d would need about 200 days to 
build up to the 0.6 cycle advance shown for 2002 in the top 



plot of Fig. 12. Then, in 2002, the peak of the periodogram 
is 1.5383 c/d, an apparently rather strong modulation and 
stable frequency, judging by the size of the error bars on 
Fig. 12 as well and the points' linearity, respectively. The 
strength can be seen in the amplitude plot. Fig. 6. 

A slight increase in frequency for both modulations oc- 
curred between the mid-90's and late-90's, as can be seen 
both in the frequency plots and the phase diagram. It is un- 
clear if a single physical mechanism modified both periods, 
but the possibility should be considered for a physical model. 
It should also be noted that the phase matched across col- 
ors. Some physical models may imply a detectable phase lag 
between band passes, but our analysis did not see any such 
lag. 



5 DISCUSSION 

iHubert Flooued ^^M) used Hipparcos photometry to dis- 
cover that short-term periodic variability (<3.5 days) was 
present in about 40% of Be stars between BO and B3 and is 
less frequent in later tvDes. llIaberl et all |2004b) calculated 
that both stars A and B ha ve spect r al typ e 09 V from the 
data in the catalog given bv lMassevI (|2002r) . so the presence 
of a short-term period is not surprising. 

Neither of these periodicities are likely to be due to 
orbits around the star. From the Corbet relation (jCorbed 
, the expected orbital period is hundreds of days be- 
cause it is correlated with the pulsar's long period of 701.6 
s. Also, if a solid body in Keplerian rotation were driving 
the variability, there would not have been glitches in phase 
stability. For the same reason neither of the variations are 
likely to be caused by rotational modulation of a surface in- 
homogeneity . although such interp retation is the object of 
controversy fBaad e L Balonal[l994r ). 

The periods should not be imputed to variations in the 
disk because the phase of the variation would not be so sta- 
ble: disks in Be stars are transient. It remains to be proven 
that there is a disk around star B; there was no coincident 
H-alpha emission in the catalog of Mevsson nier AzzoDar"dil 
(1993). The only other (simple) explanation is stellar pulsa- 
tions. 

Radial pulsations (i.e. p Cephei stars) are a possibility. 
The driving force for radial pulsation has been ident ified 
wi th an opacity (k) mechanism o f iron dCox et al.l ll992Y so 
as lSterken JerzvkiewicJ (|l988l) pointed out, it should be 
a rarity the Magellanic Clouds. I n a recent OGLE-II study 
of the Large Magellanic Cloud (jPierulski Kolaczkowskil 
2002), three P Cephei stars were found, and a study of the 
SMC for more was promised. Periods from radial pulsations 
of the J3 Cephei type in the Galaxy are generally less than 
0.3 d (jOautschv Saiolll996l ). which fits for the shorter pe- 
riod of 0.285 d but no t for t he longer period of 0.649 d. 
However, Kolaczkowski et all (|2p03) have found 20 doubly 
periodic B stars in the LMC with periods in the vicinity of 
fi and /2, so this star may be the first SMC version of those 
stars. The authors tentatively assigned those stars to both 
(3 Cephei and Slowly Pulsating B-star types. The following 
differences caution that assignment of star B to that group 
is non-trivial: this star has an earlier spectral type than any 
of those stars, it is in the SMC which has lower metallicity 
than the LMC, and it is part of a HMXB. 
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Figure 12. The data, after being cleaned of the other frequency, is spHt into chunks of about 20 observations and folded at each 
frequency. The phase with a Monte-Carlo error estimate is plotted. If the chunk has peak-to-peak amplitude less than 10 mmag, I judged 
the fit to be insignificant and did not plot a phase for that chunk. Phase difi"erence between adjacent points is forced under 0.5 cycles to 
impose continuity, but that decision is aesthetic: true phase may be an integer value ofi". Open circles are MACHO blue, filled circles are 
MACHO red, and crosses are OGLE /-band. 



© 2004 RAS, MNRAS 000,[TH8l 



8 Fabrycky 



The blue band would be most sensitive to the temper- 
ature changes of a hot, radially pulsating star, and there is 
a marked increase in amplitude at higher wavelength for /2 . 
During a cycle, the highest luminosity of such a star lags 
behind its highest temperature (£ox 1980). This effect is 
unlikely to induce a detectable phase lag between our band- 
passes, although this assertion should be quantified by future 
work. 

Many recent studies describe the short term variabil- 
ity of Be stars (and hot stars in general) in terms of non- 
radial pulsations (NRP). Such studies put an emphasis 
on time-resolved spectroscopy, looking for line-profile vari- 
ability. Photometric variations are thought to be less pro- 
nounced since the whole star and its disk co ntribute, possi- 
bly ca nceling out a locally strong variation. iR.ivinius et al.l 
(|2nn3l) give line-profile variability for 27 early-type Be stars, 
and find that / = |m| = 2 pulsations can explain nearly all 
those observations. A few stars have been monitored exten- 
sively, and multiperiodicity has been detected in a few cases. 
Ijanot-Pacheco et al.l (|l999i) found strong photometric peri- 
odicity of 1.56 c/d and intermittent periodicity of various 
frequencies, including 3.23 c/d, in 77 Cen. The k, mechanism 
has not been conclusively linked to NRP, so the low metal- 
licity of the SMC is not necessarily a problem. In line pro- 
file variations, features are seen traveling from blue to red 
(iQautschy (fe Saio|l99(t). It remains to be determined if NRP 
would cause a detectable phase lag in the our colours. The 
amplitude of variation is expected to depend on colour in 
NRP, but different parts of the star are relatively cold and 
hot simultaneously, so it is probably weaker than for radial 
pulsations. NRP seems a likely source for the /i variability 
and a possible but unlikely source for the variation at /2. 

In conclusion, I suggest that the periodic variations of 
star B are due to stellar pulsations, but which ones cause 
them is left to future work. Tentative evidence has been 
given that fi is a low-order non-radial pulsation and /2 is a 
radial pulsation. Spectroscopy and further analysis is needed 
before the cause of this multiperiodicity can be understood 
with rigor. In particular, the star's metallicity is one impor- 
tant parameter to measure, and observations of line-profile 
variability would be very valuable to discern between radial 
and non-radial pulsation. 
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